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hybrid inflation and leptogenesis, but without taking the specific details of these sce- 
narios. It is shown that the Peccei-Quinn field can successfully act as a curvaton 
generating the total curvature perturbation in the universe in accord with the cos- 



mic background explorer measurements. A crucial phenomenon, which assists us to 
achieve this, is the 'tachyonic amplification' of the perturbation acquired by this field 
during inflation if the field, in its subsequent evolution, happens to be stabilized for 
a while near a maximum of the potential. In this case, the contribution of the field 
to the total energy density is also enhanced ('tachyonic effect'), which helps too. 
The cold dark matter in the universe consists, in this model, mainly of axions which 
carry an isocurvature perturbation uncorrelated with the total curvature perturba- 
tion. There are also lightest sparticles (neutralinos) which, like the baryons, originate 
from the inflationary reheating and, thus, acquire an isocurvature perturbation fully 
correlated with the curvature perturbation. So, the overall isocurvature perturbation 
has a mixed correlation with the adiabatic one. It is shown that the presently avail- 
able bound on such an isocurvature perturbation from cosmic microwave background 
radiation and other data is satisfied. Also, the constraint on the non-Gaussianity of 
the curvature perturbation obtained from the recent Wilkinson microwave anisotropy 
probe data is fulfilled thanks to the 'tachyonic effect'. 
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1. Introduction 

The recent data |j] on the acoustic peaks of the angular power spectrum of the cosmic 
microwave background radiation (CMBR) strongly favor the idea of inflation, which 
offers the most elegant solution to the outstanding problems of standard big bang 
cosmology and exterminates unwanted relics such as monopoles. Moreover, inflation 
is now generally accepted as the most likely origin of the primordial density pertur- 
bations which are needed for explaining the structure formation in the universe [Q]. 
The usual assumption is |2|, |3|] that these perturbations come solely from the slowly 
rolling inflaton field. In this case, the observed density perturbations are expected to 
be purely adiabatic, since fluctuations in the inflaton cannot cause an isocurvature 
perturbation. Also, significant non-Gaussianity is not encountered Q in the usual 
one-field inflationary models, while the non-Gaussianity which is possible || in multi- 
field models requires extreme fine-tuning of the initial conditions. However, although 
adiabatic and Gaussian perturbations are perfectly consistent with the present data, 
appreciable non-Gaussianity || or the presence of a significant isocurvature density 
perturbation R, ffl cannot be excluded by observations. 



An alternative possibility H, which has recently attracted attention [10, 11 



[Ti2| , P"5| , I4"[ |T5"| , [TT)| , I7| , is that the adiabatic density perturbations originate from 



- 1 - 



the vacuum perturbations during inflation of some light 'curvaton' field different 
from the inflaton. In the curvaton scenario, significant non-Gaussianity may easily 
appear because the curvaton density is proportional to the square of the curvaton 
field amplitude. Also, the curvaton density perturbations can lead, after curvaton 
decay, to isocurvature perturbations in the densities of the various components of 
the cosmic fluid. In the simplest case, the residual isocurvature perturbations are 
either fully correlated or fully anti-correlated with the adiabatic density perturbation, 
with a calculable and generally significant relative magnitude. In the presence of 
axions, however, the correlation is, in general, mixed. It is important to be noted 
that the curvaton hypothesis makes |13| the task of constructing viable models of 
inflation much easier, since it liberates us from the very restrictive requirement that 
the inflaton is responsible for the curvature perturbations. 

The most compelling extension of the standard model (SM) of particle physics is 
its supersymmetric (SUSY) version, the so-called minimal supersymmetric standard 
model (MSSM). It is, however, clear that even the MSSM must be part of a larger 
scheme since it leaves a number of fundamental questions unanswered. For instance, 
the strong CP and /i problems cannot be solved in MSSM. Also, it is not easy to 
generate in MSSM the observed baryon asymmetry of the universe by electroweak 
sphaleron processes or to realize inflation. 

The strong CP problem can be elegantly resolved by including a Peccei-Quinn 
(PQ) symmetry |18|| broken spontaneously at an intermediate mass scale, which can 
be easily generated within the supergravity (SUGRA) extension of MSSM. The PQ 
field, which breaks the PQ symmetry, corresponds, in this case, to a flat direction 
in field space lifted by non-renormalizable interactions. Moreover, the // parameter 
of MSSM can be generated from the PQ scale. A minimal extension of MSSM, 
which solves the strong CP and /i problems along these lines, has been constructed in 
Ref. ||20|| . A key ingredient was a global U(l) R symmetry obeyed by the superpoten- 
tial. This model has been further extended |2l| to simple SUSY grand unified theory 



(GUT) models leading to hybrid inflation [22 and yielding successful baryogenesis 
via a primordial leptogenesis [23], [24| in accord with the data on neutrino masses and 
mixing. Also, the R symmetry implies exact baryon number conservation, thereby 
explaining proton stability. 

In this paper, we examine whether, in the above models, the PQ field, which pos- 
sesses an almost flat potential, can also play the role of the curvaton generating the 
adiabatic curvature perturbations with a possibly significant non-Gaussian compo- 
nent, while accompanied also by an isocurvature density perturbation. An important 
requirement is that the PQ field is effectively massless during inflation. In view of 
the quasi-flatness of the potential, one would imagine that this condition could be 
easily fulfilled by just ensuring that the value of the PQ field during inflation is not 
too large. However, as it is well-known [^5], large SUGRA corrections during 
inflation can destroy the quasi-flatness of the PQ potential and, thus, invalidate the 
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possibility of using the PQ field as curvaton. Therefore, it is crucial to assume that 
some mechanism [27| is employed to keep these corrections under control. 

After the end of inflation, the inflaton starts performing coherent damped os- 
cillations and eventually decays reheating the universe. We assume that the SUSY 
breaking corrections to the PQ potential which arise |28] from the finite energy 
density of the oscillating inflaton field are also relatively small. Under these circum- 
stances, the PQ field remains overdamped and, thus, undergoes a slow evolution for 
quite some time after the termination of inflation. At later cosmic times, however, 
it settles into coherent damped harmonic oscillations about zero or the PQ vacua, 
depending on its value at the end of inflation. Actually, we find that the values of 
the PQ field at the end of inflation fall into well-defined alternating bands which lead 
to the trivial (false) vacuum or the PQ vacua in turn. Of course, bands leading to 
the trivial vacuum must be excluded. 

The perturbation acquired by the (light) PQ field during inflation evolves at sub- 
sequent times and, when the system settles into damped quadratic oscillations about 
the PQ vacua, yields a stable perturbation in the energy density of the oscillating 
field. After the PQ field decays, this perturbation is transferred to the radiation 
dominated plasma, thereby generating the total curvature perturbation. We find 
that, if the PQ field, as it evolves, happens to stay near a maximum of the potential 
for some time, its primordial perturbation from inflation is drastically enhanced and, 
thus, the resulting energy density perturbation of the oscillating field can be quite 
sizeable. Also, the contribution of the oscillating field to the total energy density 
comes out, in this case, larger due to the delay in the onset of the field oscillations 
caused by the temporary rest of the field near a maximum. Thus, the final curva- 
ture perturbation can be adequate to meet the cosmic background explorer (COBE) 
results [p9[ . This 'tachyonic amplification' of the inflationary perturbation of the 
PQ field as well as the enhancement of its contribution to the energy density, which 
occur well after the end of inflation, are crucial for the viability of our scheme. 

The cold dark matter (CDM) of the universe, in our case, consists of axions and 
lightest neutralinos. The latter originate from the late decay of the gravitinos which 
are generated during the reheating process which follows inflation. In accordance to 
the curvaton hypothesis, the inflaton and, thus, the radiation which emerges from its 
decay at reheating possess practically no partial curvature perturbation. Therefore, 
the baryons, which originate from a primordial leptogenesis taking place at reheating, 
as well as the neutralinos carry negligible partial curvature perturbation. Then, at 
curvaton decay, where its curvature perturbation is transferred to the radiation, the 
baryons and the neutralinos acquire an isocurvature perturbation, which is fully 
correlated with the total curvature perturbation. The axions, which come into play 
at the QCD phase transition, also acquire [[HJ an isocurvature perturbation which 
is, though, uncorrelated with the total curvature perturbation. So, we finally obtain 
an isocurvature perturbation of mixed correlation with the curvature perturbation. 
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We apply the most stringent available restriction on such isocurvature perturbations. 
This restriction results from the set of the CMBR and other data used in Ref. |~5] . 

In our scheme, the curvaton decays somewhat before dominating the energy den- 
sity of the universe. This can generate appreciable non-Gaussianity in the curvature 
perturbation. We consider the bound on the possible non-Gaussian component of 
the curvature perturbation from the recent CMBR data which were obtained by the 
Wilkinson microwave anisotropy probe (WMAP) satellite. The 'tachyonic effect', i.e. 
the possible rest of the field near a maximum of the potential for a period of time, 
is essential also for the fulfillment of the non-Gaussianity requirement. 

The paper is organized as follows. In section |2], we summarize the minimal 
extension of MSSM which solves the strong CP and /z problems via a PQ and a R 
symmetry. We also indicate the salient features of the simple SUSY GUTs which 
further extend this model and yield hybrid inflation and successful baryogenesis. In 
section [3|, we outline the evolution of the PQ system during inflation, the subsequent 
inflaton oscillations and the period after reheating. In section [|, we discuss the 
conditions under which the PQ field can successfully act as a curvaton. In particular, 
we consider the constraints on the total curvature and isocurvature perturbation from 
the COBE results and other data. The requirements from the CDM in the universe 
are also taken into account. Section ^| is devoted to the detailed numerical analysis 
of the evolution of the PQ field and its perturbations which originate from inflation. 
The various requirements of section |] and the constraints from the non-Gaussianity 
of the curvature perturbation are numerically studied. Finally, in section we 
summarize our main conclusions. 



2. A simultaneous solution of the strong CP and /i problems 



It is well-known that, in the context of the SUGRA extension of MSSM, appropriate 
flat directions in field space can generate an intermediate scale 



Mi ~ (m 3 / 2 fflp) 



10 



ii 



10 12 GeV, 



(2.1) 



where m 3 / 2 ~ 1 TeV is the SUSY breaking scale and mp = (SttGn)^ 1 ^ 2 — 2.44 x 
10 18 GeV is the 'reduced' Planck mass with Gn being the Newton's constant. It is 
then natural to identify the scale Mj with the breaking scale f a of a PQ symmetry 
£7(1)pq [18] , which resolves the strong CP problem. Moreover, the /i parameter of 
MSSM could be obtained as \i ~ fa/ m P ~ m 3/2 [|H . 



This idea has been elegantly implemented in Ref. [2(| in the presence of a non- 
anomalous R symmetry U(l)p. To this end, an extra pair of gauge singlet chiral 
superfields N, N with the superpotential couplings 



5W 



—^—N 2 N 2 + ^-N 2 h x h 2 
2mp mp 



(2.2) 
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was added in MSSM. Here, hi and h 2 are the two Higgs doublets coupling to the 
down and up type quarks respectively, and A, j3 are dimensionless parameters, which 
can be made positive by appropriate field redefinitions. The PQ and R charges are 



PQ: iV(-l), N(l), Ml), ^(1), 
R: N{1), N(0), htiO), h 2 (0) 



(2.3) 



with the 'matter' (quark and lepton) superfields having PQ = —1/2, R = 1. The 
superpotential has R = 2. Note that global continuous symmetries such as our PQ 



and R symmetry can effectively arise [31] from the rich discrete symmetry groups 



encountered in many compactified string theories (see e.g. Ref. |3"2"|). It should be 
mentioned in passing that such discrete symmetries can also provide |J3| an alterna- 
tive solution to the fi problem. 

The part of the tree-level scalar potential which is relevant for the PQ symmetry 
breaking is derived from the superpotential term XN 2 N 2 /2mp and, after soft SUSY 
breaking mediated by minimal SUGRA, is given by |2(J 



PQ 



m 3/2 + A 



NN 



mp 



(\N\ 2 + \N\ 2 ) + I Am 3/2 \ 



N 2 N 2 



2m i 



+ h.c. 



(2.4) 



where A is the dimensionless coefficient of the soft SUSY breaking term corresponding 
to the superpotential term XN 2 N 2 /2m p. 

The potential Vj>q is minimized along the field direction where N, N have the 
same magnitude and their phases are arranged so that the last term in the right hand 
side (RHS) of Eq. (|2.4j) becomes negative. The potential then takes the form [^0 



V PQ = 2\N\ 2 ml /2 I 1 - \ A\X- 



+ A^ 



\N\* 



2m 3 / 2 mp 

For \A\ > 4, the absolute minimum of the potential is at 

\A\ + (\A\ 2 -12f - 



2 2 
m 3/2 m P / 



(2.5) 



|<jv)| = IWI 



6A 



{m 3/2 m P ) 



{m 3/2 m P ) 



(2.6) 



The vacuum expectation values (VEVs) of N, N break spontaneously the symmetry 
C/(1)pq x U(1)r down to the discrete 'matter parity' symmetry Z™ p under which the 
'matter' superfields change sign. Substituting these VEVs in Eq. ( |2.2|) , we see that 
a /i-term with \i ~ m 3 / 2 is generated. 
This scheme has been embedded 



211 in concrete SUSY GUT models which 



naturally lead to hybrid inflation |22| followed by a reheating process which satisfies 
the gravitino constraint 0]. These models also yield successful baryogenesis via a 
primordial leptogenesis |2jJ which takes place [33 non-thermally at reheating by the 
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decay of right handed neutrino superfields which emerge from the inflaton decay. 
The observed baryon asymmetry of the universe can be achieved in accord with the 
present neutrino oscillation data. 

In these models, the PQ and R charges of N, N and the MSSM superfields are 
the same as above with the exception of the PQ charge of the left (right) handed 
quark and lepton superfields which is now PQ = — 1 (PQ = 0). This fact implies 
that the discrete 'matter parity' symmetry (Z™ p ) is now not contained in C/(1)pq and 
has, thus, to be independently imposed. The VEVs of N, N break U(1)j>q x U(1)r 
completely, while Z™ p again remains exact. 

A closer look reveals that the QCD instantons break £7(1)pq explicitly to its 
Z 6 subgroup, which is then broken by (N), (N) (as in Refs. p5l |36|). This would 



lead |37| to a disastrous domain wall production if the PQ transition occurred after 
inflation. Avoidance of this disaster would then require extending the model so 
that such spontaneously broken discrete symmetries do not appear. The soft SUSY 
breaking terms break U(1)r explicitly to its Z 2 subgroup which, in the presence 
of Z™ p , is equivalent to the Z 2 under which only iV changes sign (compare with 



Refs. q(|). The breaking of this Z 2 by (N) would also be catastrophic if it took 
place after inflation. Fortunately, as we will argue in section || the N, N fields 
emerge with non-zero (homogeneous) values at the end of inflation and, thus, the 
above discrete symmetries are broken. Also, there is no PQ transition after inflation. 
Moreover, as we will explain in section |5J], the system is led to the same minimum 
of the PQ potential everywhere in space. Consequently, domain walls do not form 
and there is no need to further complicate the model. 

The fact that N, N ^ after inflation solves yet another problem. It has been 
shown that the trivial local minimum of Vpq at iV = N = is separated 

from the PQ (absolute) minimum by a sizeable potential barrier. This picture holds 
[35| , j36fl at all temperatures below the reheat temperature T reh < 10 9 GeV |34| even 
if the one-loop temperature corrections to the potential are included. Thus, after 
reheating, a PQ transition would have been impossible if the system was in the 
trivial vacuum. Such a transition could have been even more impossible during 
inflaton oscillations if the system emerged after inflation being in the trivial vacuum. 
This is due to the fact that the PQ field typically acquires [28| a positive mass 2 of 
order H 2 (H is the Hubble parameter) via its coupling to the inflaton (see section |3|). 



3. The cosmological evolution of the PQ field 

We restrict ourselves on the field direction A^* = ±iVexp [i(a + vr)/2] (a is the phase 
of A in Eq. (|2.4| )) along which Vpq is minimized and takes the form in Eq. ( |2.5| ). Ro- 
tating iV on the real axis by an appropriate R transformation, we can then define the 
canonically normalized real scalar PQ field <p — 2iV, which breaks the PQ symmetry 



- 6 - 



and generates the /x-term as discussed in section [| Eq. ( |2.5|) then becomes 



v eQ = [i - \^^; + »u^r ) ■ (») 

For \A\ > 4, this potential has a local minimum at <fi = and absolute minima at 

<0> 2 = / a 2 = ^ (|A| + vW 3 ^) m 3/2 m P (3.2) 
with f a > 0. We can shift the potential Vpq in Eq. fl3~T| ) by adding to it the constant 

V ° = lk\ + V 7 !^! 2 " 12 ) [l^l (l^l + V\A\ 2 -12) - 24] ml /2 m Pl (3.3) 

so that it vanishes at its absolute minima. 

In order to be able to study the cosmological evolution of <f), we should first 
include the SUSY breaking effects which arise [25, |2J], |2S] from the fact that the 
energy density in the early universe is finite. During inflation and the subsequent 
inflaton oscillations, SUSY breaking is transmitted to the PQ system via its coupling 
to the inflaton given by the SUGRA scalar potential. The scale of the resulting SUSY 
breaking terms is set by the Hubble parameter H and, thus, these terms dominate 
over the terms from hidden sector SUSY breaking as long as H > m 3 / 2 - 

We should point out that, during the era of inflaton oscillations, the PQ potential 
acquires temperature corrections too. They originate from the 'new' radiation []39| 
which emerges from the decaying inflaton and has temperature T ~ {HmpT^ 



rch/ 



a/4 



Their scale is set by the thermally induced mass for the PQ field. A typical con- 
tribution (~ f3NT/rrip) to this mass is given by the second term in the RHS of 
Eq. ( |2.2| ) with the Higgsino lines connected to a loop via a mass insertion. However, 
this contribution exists only under the condition that the effective Higgsino mass 
(~ /3N 2 /mp) does not exceed T. For our input numbers (see section ^ and almost 
until the end of the period when H > 7713/2, we can show that, whenever this con- 
dition is satisfied, we also have (3NT/mp < H, which means that the temperature 
corrections are overshadowed by the SUGRA ones. At later times, N becomes any- 
way adequately small and, thus, the thermally induced mass for the PQ field comes 
out much smaller than H. We, thus, conclude that the temperature corrections can 
be neglected throughout the era of inflaton oscillations. 

After reheating, where H <ti m 3 / 2 as a consequence of the gravitino constraint 
(^reh i$ 10 9 GeV) [Q, the (one-loop) temperature corrections to the PQ potential 



could be the main source of SUSY breaking. However, as shown in Refs. [po[ , |3qj , 
these corrections are less important than the ones from the hidden sector. This is un- 
derstood from the fact that N, N possess only non-renormalizable interactions which 
are suppressed by mp and carry relatively small dimensionless coupling constants. 
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In our case, the SUSY breaking corrections to the PQ potential during the in- 
flationary and oscillatory phases can generally be [28| of two types: 



H 2 m P f^/m P ), (3.4) 

Hm P g(<f) 4 /m P ), (3.5) 

where /, g are some functions depending on the model, the Kahler potential and the 
phase. To the leading approximation, the correction in Eq. ( |3.4j ) yields a mass term 
<5m 2 2 /2 for the curvaton 0. For a general Kahler potential, Sm^ ~ H 2 with either 
sign possible. However, for specific (no-scale like) forms of the Kahler potential, this 
induced mass 2 might be (partially) cancelled |27| . For simplicity, we will approximate 
the correction in Eq. ( |3.4|) by the mass term. We will also assume that Srn 2 ^ is positive 
and somewhat suppressed. We, thus, write 

5rri0 = 7#, (3.6) 

where 7 ~ 0.1 and can have different values during inflation and field oscillations. 
The correction in Eq. (|3.5| ) is of the A term type with A ~ H and, although it is 
generally present, we will not take it into account in order to simplify the analysis. 
We will, thus, consider the following full scalar potential 

V{(j>) = V PQ + ^ 2 H 2 (P 2 + V . (3.7) 

The evolution of the field is governed by the classical equation of motion 

+ 3F0 + V\<j>) = 0, (3.8) 

where overdots and primes denote derivation with respect to the cosmic time t and the 
field respectively. Here, the damping term, proportional to the Hubble parameter, 
arises from the expansion of the universe. For large values of 0, the potential Vpq in 
Eq. ( |3.ip is dominated by the non-renormalizable 6 -term. Moreover, for H ^> 777.3/2, 
the induced mass in Eq. ( |3.6|) can also be important. Hence, if happens to be large 
during some stage of inflation, we can write 

V{4>) * \l 2 Hl^ + A 2 ^ + V , (3.9) 

where H in Q is the inflationary Hubble parameter. It is natural to assume that, at the 
onset of inflation, V((p) < V^, where Vj n fl is the (almost) constant vacuum energy 
density which is provided by the inflaton and drives the exponential expansion of the 
universe. Thus, <ft ma Y initially be quite large. For such values of 0, V" is dominated 
by the 6 -term in the potential in Eq. ( |3.9| ) and the system is initially underdamped 
{V" ^> H 2 n{i ). Consequently, decreases rapidly and V" quickly falls below H 2 na 
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with the system entering into a regime of slow-roll evolution. Soon after this, the 
induced mass term in Eq. (|3.9|) becomes dominant. 

The classical equation of motion holds until the quantum perturbations of 
generated by inflation, 5<p = H in f[/2TT, overshadow its classical kinetic energy den- 
sity Pkin = 2 /2. These perturbations contribute a 'quantum' kinetic energy den- 
sity of order (5<fi/5t) 2 /2, where 5t ~ H.~ ft . Therefore, the slow-roll regime will end 
when pkin ~ From the slow-roll equation ~ —V'{(j))/?>H- Ul ^ derived from 

Eq. (|3.8|) , we then find that this happens |4(J at ~ 0q, where 0q is given by 

3/7 3 

H0q) * 7 2 ^i 2 nfl 0Q ~ -2^. (3.10) 

The subsequent evolution of is controlled by its quantum perturbations which 
'kick' it around to perform a 'random walk'. This results to a 'quantum drift' of 
the average value of towards smaller values (since it cannot be pushed over 0q) 
in the sense that the field's distribution spreads towards the origin, which is finally 
engulfed inside an evenly spread stationary distribution whose width is 0q. Thus, 
the value 0/ of at the end of inflation, taken positive without loss of generality, 
can be smaller than 0q if sufficiently many e-foldings take place after 0q, but it is 
typically ~ 0q. At the termination of inflation, the quantum perturbations of 
cease and the field emerges with vanishing kinetic energy. 

Of course, depending on the total number of e-foldings during inflation and the 
value of at the onset of inflation, the inflationary phase may be terminated while 
is still in the slow-roll regime. In this case, we have 0/ > 0q and 0/ — — 7 2 -ff; n fl0//3 — 
\ 2 4> 5 f/16m 2 pH in{ [. For the values used here, the corresponding kinetic energy density 
is much smaller than the potential energy density and can be neglected. In summary, 
the field can practically have any value 0/ at the end of inflation. However, if 0/ 
is smaller than a certain large limit (depending on A) so that is slowly rolling or in 
the 'quantum' regime as inflation ends, the field emerges with vanishing velocity. 
This sets the initial conditions for the subsequent evolution of the PQ field. 

After inflation, the universe enters into a matter dominated era where the energy 
density is provided by the coherent oscillations of the inflaton. The potential V(0) is 
then given by Eq. fl3.7| ) with H = 2/3t. If is slowly rolling at the end of inflation, 
the friction term in Eq. ( |3.8|) remains important at subsequent times too and the 
evolution of is slow for quite some time. As the cosmic time increases, the Hubble 
parameter decreases and, thus, at some point, the induced mass in Eq. ( |3.6|) becomes 
< m.3/2. After this point, the soft terms arising from SUSY breaking by the energy 
density of the oscillating inflaton are smaller than the ones from hidden sector SUSY 
breaking. Therefore, the potential can be approximated as 

V{4>) — Vpq + V , (3.11) 
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with Vpq given in Eq. Q3.1J) , and has a local minimum at the origin with V = Vq, 
absolute minima at = ±/ a with V = and local maxima at = ±0 max , where 

4L* = ^ (\ A \ - V\A\ 2 -12) m 3/2 m P (3.12) 

with 

Knax = ^(l^l 2 " l2)lml /2 m P (3.13) 

and max > 0. Soon after this, H becomes < m 3 / 2 and the system enters into a phase 
of damped oscillations about either the trivial minimum at = or one of the PQ 
minima at = ±/a, depending on the value of at t ~ ^3/2 • 

After reheating, if = 1 / 2t -C m 3 / 2 . So, the system remains underdamped and 
V{4>) is still given by Eq. ( |3.11| ). (Recall that the temperature corrections to the 
potential are not important.) The field keeps performing damped oscillations 
about one of the minima until it finally decays. These minima are separated by 
sizeable potential barriers and, thus, a transition from the trivial to a PQ minimum 
cannot occur. So, if the field was driven to the trivial minimum, it would be trapped 
there forever. The PQ symmetry would then remain unbroken and the strong CP 
problem would not be solved, while the universe would experience an unacceptably 
large effective cosmo logical constant. 



4. The PQ field as curvaton 



In accordance with the recently investigated jR], [LI], [R| [R| 0, 0] curvaton 
scenario, the curvature perturbations may originate from the quantum perturbations 
during inflation not of the inflaton but rather of another light field which has been 
called curvaton. In this section, we will address the question whether our PQ field 
can successfully act as a curvaton. To this end, we must first consider the following 
two requirements: 

1. As already mentioned, the energy density of the universe during inflation must 
be dominated by Vj n fl, the vacuum energy density provided by the inflaton. 
(Here, V- in ^ is assumed to be essentially constant.) Therefore, the value of the 
PQ field during and at the end of inflation should be bounded above by the 
requirement V((f>) < V in R. For definiteness, we take here the standard SUSY 
hybrid inflationary scheme |?6|, [H[ (for a recent review see e.g. Ref. []36]]). The 
inflationary scale is then given by 

V& = i£m, (4.1) 

where k and M are the dimensionless parameter and the mass scale of the 
standard superpotential for hybrid inflation respectively. For simplicity, we 
identify M with the SUSY GUT scale M GUT ~ 2.86 x 10 16 GeV. 
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2. A much more stringent constraint is obtained from the requirement that the 
field is effectively massless during (at least) the last 50 — 60 inflationary 
e-foldings so that it receives a superhorizon spectrum of perturbations |42[ 



5 * = If ^ 

from inflation. This requirement means that V"((j)) < i J? fl , which guarantees 
also that is of negligible time derivative at the end of inflation. 

The detailed evolution of the inflationary perturbations of after inflation will 
be considered in the next section. For the moment, we will only discuss how these 
perturbations can generate (or affect) the total curvature perturbation 

C = -H*B- = , 6p (4.3) 
P 3(p + p) 

where p and p are the total energy density and pressure in the universe respectively. 
As we approach the time at which the curvaton decays, the universe is normally 
dominated by radiation and the oscillating curvaton 0. In this case, we can write 

C = (1 - fXr + fa, (4-4) 

where ( r = 5p r /4p r , Q = Sp^/Sp^ and 

3p4, + 4p r 

with p r and p^ being the radiation and energy densities respectively. Here, we 
assumed that the amplitude of the oscillating has been adequately reduced so that 
the potential can be approximately considered as quadratic. The oscillating curvaton 
then resembles pressureless matter and Sp^/p^ = 2<50 O /0 O , where O is the amplitude 
of the oscillations and 50o the perturbation in this amplitude originating from <50 
in Eq. (|4.2|) . According to the curvaton hypothesis, the inflaton contribution to the 
density perturbation is insignificant. Consequently, the radiation emerging from the 
inflaton decay has Q r ~ 0. After the curvaton decay, we obtain the final curvature 
perturbation which, in the instantaneous decay approximation, is given by 

C *L ideeC/,, (4.6) 

where /dec is the value of / at the decay of the curvaton. The COBE measurements 
require p9j that ( ~ 5 x 10~ 5 . 



If the curvaton field decays well before dominating the energy density of the 
universe (i.e. /dec "C 1), which is the case here, the total curvature perturbation can 
acquire [DJ a significant non- Gaussian component. The reason is that the curva- 
ture perturbation is proportional to the perturbation of the curvaton energy density 
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(see Eq. (L6)), which depends on the square of (p. This is interesting since the 
recent CMBR data from the WMAP satellite leave || considerable room for non- 
Gaussianity in the curvature perturbations. We will come back to this point in the 
next section where a detailed numerical analysis of the model will be presented. 

As mentioned in section |3], after the end of inflation, the inflaton performs 
damped oscillations and eventually decays into light particles reheating the universe. 
The reheat temperature T re h is related || to the inflaton decay width r in fl by 

/ 45 \ * i 
T reh = U-2— (r infl mp)5, (4.7) 
yZ7T g* J 

where g* is the effective number of degrees of freedom. After reheating, gravitinos 
(g) are thermally produced besides other particles. They only decay, though, well 
after the big bang nucleosynthesis (BBN) since they have very weak couplings. The 
decay channels normally include the one to a photon and a photino (g — > 7 + 7) with 
the emerging high energy photons causing the dissociation of the light elements and, 
thus, destroying the successful predictions of BBN. This yields |[4| the strongest 
possible upper bound on the reheat temperature, which is T re h < 10 9 GeV if the 
branching ratio to photons equals unity and the gravitino mass m 3 / 2 is in the range 
of several hundreds of GeV. However, the upper bound can be considerably relaxed 
[34] , |4"3|| if this branching ratio is less than unity. 

Each decaying gravitino yields one sparticle subsequently turning into the lightest 
sparticle (LSP), which is stable. This is due to the presence of 'matter parity' in 
our model. Moreover, no pair annihilation of the resulting LSPs is possible since the 
universe is already quite cold. Consequently, these LSPs survive until the present 
time contributing to the relic abundance of CDM in the universe. We will assume, 
as it is customary, that the LSP is the lightest neutralino (x). We will also make the 
simplifying assumption that the relic abundance of the thermally produced LSPs is 
negligible, which holds in many cases. So, the relic abundance, Q^sph 2 , of the LSPs 
comes solely from the decaying gravitinos and can be estimated following Ref. [f43||. 
Here, ^lsp = Plsp/Vc with plsp (Pc) being the LSP (critical) energy density and 
h ~ 0.71 is the present Hubble parameter in units of 100 km sec -1 Mpc -1 . 

Our model contains axions which can also contribute to the CDM of the universe. 
Their relic abundance, £l a h 2 , can be calculated by applying the formulae of Ref. [H3 



(Q a = Pal Pc with p a being the axion energy density). For simplicity, we take the 
QCD scale A = 200 GeV and ignore the uncertainties in these formulae. Also, the 
anharmonic effects are negligible in our case and the possibility of entropy production 
after the QCD transition is not considered. The total CDM abundance is the sum 
of the relic abundance of the LSPs from the gravitino decays and the relic axion 
abundance (^cdm = Pcdm/Pc with Pcdm being the CDM energy density): 

1.175 



^LSP 



tW* 0.0074 (^)( I J^_) + ^( I5; 4^ 1 , (4 
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where ttt-lsp is the LSP mass and 6 the so-called initial misalignment angle, i.e. the 
phase of the complex PQ field at the end of inflation. The angle 9 lies in the interval 
[— 7r/jV, ir/Af], where J\f is the absolute value of the sum of the PQ charges of all 
fermionic color (anti)triplets. This same Af determines the subgroup of [7(1)pq 
which remains unbroken by instantons. As mentioned in section |2|, in our (and the 
simplest) model, M = 6. All #'s in the above interval are equally probable. 

It is well-known that all the effectively massless (non conformally invariant) 
fields acquire a superhorizon spectrum of perturbations during inflation. So, such 
perturbations are generated in the axion field too. When the axion mass comes into 
play at the QCD phase transition, these perturbations become isocurvature axion 
density perturbations. The late universe consists, in our case, of axions (a), baryons 
(B) and neutralinos (x) as matter, and photons (7) and (massless) neutrinos {u) 
as radiation. In order to extract the restrictions on our model which come from 
isocurvature perturbations, we employ the analysis of Ref. |§. We, thus, consider 
the comoving curvature perturbation 7£ ra d and the isocurvature perturbation S ra d 
determining the large scale CMBR temperature perturbation. They are given by 

ftrad = -C, (4-9) 

5 rad a = T^-SiOi, (4.10) 

*=a,B,x 

where Si = Srii/rii — 3 5T/T (i = a, B, x) is the isocurvature perturbation in 
the ith species, a, a, are appropriate Gaussian random variables, Qi = Pi/p c and 

= Pm/Pc with pi and Hi being the energy and number densities of the ith species 
and p m = p a + Pb + P\ being the matter energy density. 

The neutralinos appear as decay products of the gravitinos which were gen- 
erated at reheating. They, thus, inherit the same curvature perturbation, which 
coincides with that of radiation emerging from the inflaton decay. However, ( r ~ 
by the curvaton hypothesis and, consequently, ( x = 5p x /?>p x also vanishes. After 
the curvaton decay, the neutralino isocurvature perturbation S x = 3(( x — ( y ) ~ —3(. 
Here, we assumed that the neutrino isocurvature perturbation is negligible and, thus, 
C 7 = 5p 7 /4p 7 ~ ( (p 7 is the photon energy density). Baryons can be produced via a 
primordial leptogenesis |23| which can occur |24[] at reheating. The inflaton decays 



into right handed neutrino superfields which, subsequently, decay into Higgs and lep- 
ton superfields generating a net lepton number. As a consequence, £g — too and 
again Sb = 3(Cb — Cy) — — 30 We see that the baryons and the neutralinos carry an 
isocurvature perturbation which is fully correlated with the isentropic perturbation 
in the notation of Ref. ||. 

The axions are produced at the QCD phase transition which takes place com- 
fortably after the curvaton decay. At this transition, the axion mass m a is gener- 
ated and the axion field starts oscillating. For the amplitudes a encountered here, 
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the anharmonic effects |PJ] are negligible and, thus, the axion energy density is 
given by p a = m 2 a a 2 /2. The isocurvature perturbation which is generated is then 
S a = 5n a /n a = 5 p a / p a = 2 5a/ a. The perturbation 5a in the amplitude of the axion 
field comes from inflation and can be estimated jEJ as follows. Suppose that 0* is the 
value of when our present horizon scale crosses outside the inflationary horizon (i.e. 
~ 50—60 e- foldings before the end of inflation) and 9 is the initial misalignment angle, 
i.e. the phase of the complex PQ field during (and at the end of) inflation. This angle 
corresponds to the canonically normalized real scalar field 9(f)* which acquires a per- 
turbation if in fl/27r from inflation and, thus, 59 = H-^/I-k^*. The initial amplitude 
of axion oscillations a = 9f a then acquires a perturbation 5a = H infl fa/2TT(f) !l: . We 
see that the axion isocurvature perturbation is S a = Hi n ^/^9(j)*. This perturbation 
is completely uncorrelated with the curvature perturbation. 
Replacing S a , Sb and S% in Eq. ( |4.10| ), we obtain 



<Sretd a 



Tit 



a B 



(4.11) 



which consists of an uncorrelated and a fully correlated part. Here, the Gaussian 
random variables a a , a# = a^ are independent. The dimensionless cross correlation 
parameter is given by H 



cos A = — 
B 



Vt m h 



(4.12) 



where 



B 



infl 



Tit 



+ 9 1 



ttgh 2 



(4.13) 



is the entropy-to-adiabatic ratio. In Ref. ||, a bound on this ratio has been derived 
from the pre-WMAP CMBR data and was depicted in Fig. 2 of this reference versus 
cos A with the other parameters marginalized. We take an improved version |4lJ of 
this figure, which can be derived from the extended set of pre-WMAP data used in 
Ref. |p.5|| , and extract from it the 95% confidence level (c.l.) bound on B for each 
value of cos A. At present, this is the most stringent available bound on isocurvature 
perturbations of mixed correlation with the curvature perturbation. 



5. Numerical analysis and results 
5.1 The band structure 

We are now ready to study numerically the evolution of the PQ field after inflation. 
To this end, we perform the change of variable 

t = Hry (5.1) 
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in Eq. (|3.8|) , which then becomes 

d 2 <f) , f o T H_ _ 1 \d0 + e? T 
hn& J dr Hr Qii 



Inflation ends at t = tf ~ (2/3)^"^. So, we impose the initial conditions = 0/ > 
and d(f)/dr = at r = ln(2/3). Although, as explained in section [| the field exits 
inflation with negligible velocity only if 0/ is not too large, we will apply these initial 
conditions uniformly for all values of 0/ in order to obtain an overall picture of the 
possible behavior patterns of the system. We will, of course, use only the relevant 
0/'s to draw physical conclusions. As we already mentioned, these 0/'s practically 
coincide with the 0/s for which the PQ field is effectively massless during inflation 
so that it acquires a superhorizon spectrum of perturbations. 

The Hubble parameter during inflation H iri Q is given by the Friedmann equation 



n 2 _ ^nfl _ k M GUT _ o 

3mp 3mp 



where the almost constant energy density driving inflation is taken from Eq. (|4T 
which holds in the case of standard SUSY hybrid inflation. We do not commit 
ourselves to any details of this inflationary scheme. We also leave out the possible 
restrictions from the specific reheating and leptogenesis processes in the SUSY GUT 



models of Ref. pij . We just use, for our analysis, the inflationary scale derived from 
Eq. ( flip as a working example. For definiteness, we further take k = 1(U 3 , which 
gives V^ 4 ~ 9.04 x 10 14 GeV and H mR ~ 1.97 x 10 11 GeV. It should be noted that 
this inflationary scale is consistent with the curvaton hypothesis which requires [ITS 



that V^Jfi < 3 x 10 15 GeV. Indeed, if this condition holds, the curvature perturbation 
from the inflaton can be made negligible. 

The scalar potential V(0) is taken from Eq. (|3.7|) , where we put m 3 / 2 = 300 GeV, 
\A\ = 5 and 7 = 0.1. The parameter A is taken to range between 10~ 4 and 10 -3 . 
Fixing the reheat temperature T reh to the value 10 9 GeV |54j and taking = 228.75 
which corresponds to the MSSM spectrum, Eq. (|4.7|) yields the inflaton decay time 
£ re h = ~ 2.44 x 10 _1 GeV -1 . For tf < t < t re h, the universe is dominated by the 
oscillating inflaton and H = 2/31 The expression in parentheses in Eq. ( p.2[ ) is then 
equal to unity. For t > t reh , however, radiation dominates and H = l/2(t — t re h/4). 
Here, we subtracted leh/4 from t in order to achieve continuity at t = t Te ^. For 
t ^> t re h, the expression in parentheses in Eq. (|5~2] ) equals 1/2. 

Using the above described initial conditions and numerical inputs, we now solve 
Eq. ( |5.2|) numerically. For each value of <pf and any given A in the range 10~ 4 — 10~ 3 , 
we follow the evolution of the system and identify the minimum of Vpq about which 
the PQ field is finally set in coherent damped oscillations. We find that a well- 
defined band structure emerges in the A — 0/ plane. This structure is depicted in 
Fig. [I| with alternating white and green bands. We number the successive bands 
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Figure 1: The band structure in the A — (j)f plane. If 4>f lies in a green (white) band, the 
system is led into a PQ (the false) vacuum. The red crossed area is excluded by requiring 
that the PQ field is light during the relevant part of inflation, while the red dotted area is 
excluded by the bound on the isocurvature perturbation from the CMBR data. 



starting from zero which is given to the lowest lying (white) band. So, the white 
bands correspond to even integers, while the green ones to odd integers. We only 
show the first 18 bands since after this the bands become too dense in the figure to 
be clearly distinguishable. If the initial value of the field <pf belongs to a white band, 
the system is finally led into the false vacuum at cf) = 0. So, such <^'s are excluded. 
However, if 0/ lies in a green band, the system ends up in one of the PQ vacua. The 
boundaries of the bands in Fig. BL which is a log — log plot, are parallel straight lines 
with an inclination equal to —1/2. This is easily understood by observing that, if <fi 
is rescaled by absorbing a A 1//2 factor, Eq. (|5.2j) becomes A-independent. 

As indicated in section ||, the PQ field must be light during the last 50 — 60 
e-foldings (i.e. V"(</>) < #? fl ) so that it acquires perturbations from inflation which 
enable it to act as curvaton. In particular, at the end of inflation, we must require 
that V"(<pf) < Hi n R, which excludes the red crossed area in Fig. |I[ The boundary of 
this area is again a straight line with an inclination equal to —1/2 and turns out to 
be very close to the lower boundary of the 5th (green) band. So, out of all the odd 
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(green) bands, only the lrst and 3rd are consistent with the above requirement, while 
all the higher green bands are excluded. We find that, in the non-excluded area, 
remains practically constant in the last 50 — 60 e-foldings and, thus, the condition 
V'(</>) < H? nfl is valid throughout this period. It is obvious that, in this same area, 
the field velocity at the end of inflation is negligible. Finally, we should note that, in 
this allowed area, the condition V{<p) < V- m ^ which, as explained in section f|, must 
hold during inflation, is automatically satisfied provided that the total number of 
inflationary e-foldings is not huge. 

It should be emphasized that the non-vanishing of N, N after inflation alone 
cannot guarantee the avoidance of disastrous domain walls. It is crucial to further 
ensure that the system finally relaxes in the same minimum over all space. This is 
readily fulfilled in our case since the perturbation 6(f) = H inR /27r ~ 3.14 x 10 10 GeV, 
which is acquired by the PQ field during inflation, is much smaller than the width 
of the bands. Therefore, right after inflation, the field is found to lie within the 
same band everywhere in space. It is then led into the same minimum over all space 
and no catastrophic domain wall production occurs. The opposite situation has 
been encountered in Ref. fl47 |) where the inflationary perturbation of the field could 
become comparable to or even larger than the width of the bands. As a consequence, 
topological defects could be generated. 

We will now consider the bound on the entropy-to-adiabatic ratio B in Eq. ( |4.13j ) 
which has been derived in Ref. || from the pre-WMAP CMBR data. We take the 
baryon abundance QbIi 2 = 0.02, which is |48| its central value from BBN. The CDM 
abundance, f2 CDM /i 2 , is taken equal to its central value 0.12 from DASI JEJ. Thus, 
the total matter abundance is Q m h 2 = 0.14. The axion abundance fl a h 2 (second 
term in the RHS of Eq. (|4.8| )) can then be calculated from Eq. fl4.8|) , where we put 
^lsp — 200 GeV. We find that the axions constitute about 94% of the CDM. From 
this, we find the initial misalignment angle 9 for each A. The curvature perturbation 
( is taken equal to 5 x 10~ 5 from COBE p9| . As we saw, <fi remains practically 
constant in the last 50 — 60 e-foldings provided that we are outside the red crossed 
area of Fig. [l|. Thus, in the relevant area, the value 0* of <fi when our present 
horizon crossed outside the inflationary horizon can be identified with <fif. We apply 
an improved version |46| of the 95% c.l. bound on B for each value of the cross 
correlation parameter cos A in Eq. ( f4.12| ). This yields a lower bound on 0j for each 
A excluding the red dotted area in Fig. [TJ. The boundary of this area is a straight 
line with an inclination equal to —1.175/4 ~ —0.294 as one can easily deduce from 
Eqs. (pD, (p|), and (fO|. 

We see that the bound from the isocurvature perturbation excludes the lrst 
(green) band and, thus, only the 3rd (green) band survives. Note, though, that 
this bound actually pushes 0* and not <pf to large values. Moreover, in the last 
50 — 60 e-foldings, such a large 0* could be reduced to a <pf lying in the lrst band 
if <p was rolling fast. However, this would require a large SUGRA induced curvaton 
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mass during inflation, which is incompatible with the curvaton hypothesis. Finally, 
we should stress that the fact that the CDM consists mainly of axions carrying an 
isocurvature perturbation which is uncorrected with the curvature perturbation is 
of crucial importance for the viability of the model. Indeed, it leads to a cross 
correlation parameter cos A which is somewhat smaller than unity and, thus, the 
upper bound on the entropy-to-adiabatic ratio B is considerably relaxed. Actually, 
in our case, this bound takes its maximal value (~ 0.95) which is achieved in the 
range 0.55 < cos A < 0.75. 

5.2 The evolution patterns of the PQ field 

Let us now describe in some detail the time dependence of the full curvaton potential 
V((p) in Eq. ( |3.7p and the evolution of the PQ field after the end of inflation in 
accordance with our numerical findings. For definiteness, we fix A = 10 -4 . To obtain 
the picture for other values of A, we just have to rescale V as A -1 and as A -1 / 2 . 
In Fig. we present a three dimensional plot of V{<p) as a function of and the 
cosmic time t. We see that, initially and for small 0's, the induced mass in Eq. ( |3.6| ) 
dominates over the soft SUSY breaking terms in Eq. (|2.4[) and, thus, V((p) has just 
one minimum at zero. At later times, the induced mass becomes comparable to 
the soft terms and V(4>) develops a pair of symmetric local minima separated from 
the trivial minimum (at zero) by a pair of symmetric local maxima. The height 
of the local minima is gradually reduced and, at some point, they become absolute 
minima. Finally, when t > 2-f/3m 3/2 ~ 2.22xl0~ 4 GeV -1 , the induced mass becomes 
subdominant and, thus, V(4>) can be approximated as in Eq. ( |3.11| ). The above 



(non-trivial) minima then approach the PQ vacua at = ±f a ~ ±6.48 x 10 12 GeV 
with V = 0. The height of the local maxima also decreases as the universe evolves 
and, for t > 27/3m 3 / 2 , they approach the maxima of the (approximate) potential 
in Eq. fl3.11| ), which lie at = ±0 ma x ~ ±2.61 x 10 12 GeV with potential energy 
density V = V max ~ (3.27 x 10 7 GeV) 4 (see Eqs. (gT2l) and (gTgD). The trivial 



minimum (at = 0) exists at all times and possesses a constant potential energy 
density V = V ~ (3.16 x 10 7 GeV) 4 (see Eq. (p|)). 

In order to get a feeling of the post-inflationary evolution patterns of in the var- 
ious bands, let us first consider that 0/ is slightly above (below) the lower boundary 
of the lrst band which lies at 0/ ~ 2.83 x 10 12 GeV for A = 10 -4 . The inter- 
esting part of the path followed by is depicted in Fig. ga by a green (red) line. 
The field is initially overdamped and decreases monotonically at a very slow pace. 
For t > 27/3m 3 / 2 ~ 2.22 x 10 -4 GeV -1 , is found to lie slightly above (below) 
max ~ 2.61 x 10 12 GeV. For t > 2/3m 3/2 ~ 2.22 x 10~ 3 GeV -1 , H < m 3/2 
and the system becomes underdamped. It falls towards the PQ vacuum at = f a 
(the trivial false vacuum) and starts performing damped oscillations about it. Soon 
after this, the universe enters into the radiation dominated era which follows reheat- 
ing and the damped oscillations of continue until this field finally decays. For 
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Figure 2: The three dimensional plot of the full curvaton potential V{4>) versus 4> and t for 
A = 10~ 4 . The green (red) line in (a) depicts the path followed by (ft if 4>f is a little above 
(below) the lower boundary of the lrst band, while the green (red) line in (b) corresponds 
to the field evolution if 4>f is a little below (above) the upper boundary of the same band. 

cf)f < 2.83 x 10 12 GeV, which corresponds to the zeroth (white) band, the field al- 
ways rolls monotonically towards the false vacuum at zero and, eventually, is set in 
damped oscillations about it. 

In Fig. ^b, we show the red (green) path followed by the field 0, if 0/ is taken to lie 
slightly above (below) the upper boundary of the lrst band at (fif ~ 3.11 x 10 13 GeV 
for A = 1CT 4 . We see that the field, as it decreases slowly, passes from the PQ 
vacuum at <f> — f a and then climbs up the potential until it reaches a value slightly 
below (above) max - It stays there for some time and, finally, falls towards the trivial 
vacuum (the PQ vacuum at <p = f a ) and performs damped oscillations about it. For 
2.83 x 10 12 GeV < (f)f < 3.11 x 10 13 GeV, which corresponds to the lrst (green) band, 
the field always rolls slowly directly into the PQ vacuum at (j) — f a and oscillates 
about it. We should note that the upper limit of the 'quantum' regime <pQ, which 
turns out to be ~ 9.41 x 10 12 GeV (see Eq. (|3.1(J| )), lies is this band for the range of 
A considered here. 

In Fig. 0, we present the evolution of the field when 0j happens to lie near 
the lower boundary of the 3rd (green) band at <pf ~ 2.5 x 10 16 GeV again for 
A = 10 -4 . The green (red) line corresponds to 0/ being slightly above (below) 
this boundary. The earlier evolution of the PQ field is depicted in Fig. |3|a. We 
see that, well before the appearance of the non-trivial minima, the field rolls down 
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Figure 3: The PQ field evolution for <f>t near the lower boundary of the 3rd (green) band 
and with A = 10~ 4 . The green (red) line corresponds to <j)f being slightly above (below) 
this boundary. The earlier evolution is shown in (a), while the later one in (b). 



slowly, reaches the trivial vacuum and then passes to negative values slowly climbing 
up the potential. It is, eventually, stabilized for some time slightly after (before) 
— 0max and finally falls into the PQ vacuum at <fi = — f a (trivial vacuum) as is 
obvious from Fig. |]b, where the subsequent evolution is shown. We find that, for 
3.11 x 10 13 GeV < (f)f < 2.5 x 10 16 GeV, which corresponds to the 2nd (white) band, 
the field always falls into the trivial vacuum and oscillates about it. 

Fig. |] shows the green (red) trajectory followed by for (fif slightly below (above) 
the upper boundary of the 3rd band at <pf ~ 4.48 x 10 16 GeV for A = 10~ 4 . We see 
that the slowly moving field passes from zero, overtakes —f a , falls back to the PQ 
minimum at = —f a , which has appeared by that time, climbs up to a point slightly 
below (above) — max and stays there for a period of time. It then falls towards the 
PQ minimum at <fi = — f a (the trivial minimum) and is set in damped oscillations. 
Our numerical results show that, for 2.5 x 10 16 GeV < 0/ < 4.48 x 10 16 GeV, which 
corresponds to the 3rd (green) band, the field always falls into the PQ vacuum at 
<t> = —f a and oscillates about it. More complicated patterns of evolution of the PQ 
field are encountered in the higher bands, but we will not discuss them in any detail 
since these bands are both physically meaningless and useless for our purposes as 
already explained. 

It is interesting to note that the 2nd (white) band extends over almost three 
orders of magnitude in (fif and, thus, pushes the 3rd (green) band to large values of 
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Figure 4: The PQ field evolution for 0y near the upper boundary of the 3rd (green) band 
and with A = 10 -4 . The green (red) line corresponds to <j)f being slightly below (above) 
this boundary. The earlier evolution is shown in (a), while the later one in (b). 



<pf. Also, the lower bound on <fif from the isocurvature perturbation is quite sizeable. 
The combined effect of these two facts is that the 0/s which can be useful for our 
purposes are fairly large. As a consequence, the relative perturbation, 50/0, acquired 
by the curvaton during inflation is somewhat small. This leads to a small partial 
curvature perturbation, carried by the oscillating PQ field at the time of its decay. 
So, as can be seen from Eq. (|4.6|) , a relatively large /dec is required in order to achieve 
the COBE constraint on the total curvature perturbation, which makes our task a 
little tougher. The situation would have certainly been less tight if the lrst (green) 
band could be used. 

One way to understand the fact that the 2nd band comes out so wide is the 
following. Consider the classical equation of motion in Eq. ( PTS| ) during the last part 
of the period of inflaton oscillations where the induced mass in Eq. (|3.6|) can be 
neglected. For small enough 0's, the potential V(<ft) can be approximately taken to 
be quadratic and this equation becomes of the Emden type 

0+^0 + m2 /2 = O. (5.4) 

The solution of this equation is 



= — smm 3 / 2 (t - t ), 



(5.5) 



where A and to are constants depending on the initial conditions. For A < and 
< m 3 / 2 to ~ 1, the field initially (i.e. for t < m^} 2 ~ to) decreases approximately 
as — y4sinm3/2to/t. However, as t exceeds m~^} 2 ~ to, it is set in damped oscillations 
about zero. On the other hand, for A < and < m 3 / 2 to <C 1, the field initially 
decreases and reaches zero at t = t . It then passes to negative values and is stabilized 
for a while (i.e. for t < t < aroun d Am 3 / 2 < 0. It, finally, falls back into 

zero and starts performing damped oscillations about it. The former behavior is 
encountered when lies in the lower part of the 2nd band, while the latter when <pf 
reaches the upper part of this band. We see that, in the latter case, the field 0, as 
it climbs up the potential towards the maximum at = — max , 'breaks' and enters 
for a while in a plateau. This makes it more difficult for to reach this maximum. 
Thus, the upper boundary of the 2nd band is pushed to higher values. 

5.3 The inflationary perturbation of the curvaton 

We are now ready to discuss the evolution of the inflationary perturbation of 
and estimate numerically the resulting partial curvature perturbation ^ at curvaton 
decay. For any value of A in the range 10~ 4 — 10~ 3 , we take two values of at the end 
of inflation: 0j and <pf + 5(f) (5<p = both lying in the 3rd band and follow 

numerically their evolution. (Recall that <50 is tiny compared with the width of this 
band.) For both these initial values, the field ends up performing damped oscillations 
about the PQ vacuum at = — f a . When the amplitude 0o of these oscillations is 
adequately reduced, the oscillations become harmonic with their amplitudes differing 
by an amount S(J)q. The partial curvature perturbation is then stabilized to the 
limiting value 2£0o/30o, as shown in section This occurs well before the curvaton 
decay in all the cases which we considered. We find numerically the limiting for 
10~ 4 < A < and any <pf in the 3rd band. 

It should be mentioned in passing that, for a scalar field with a positive power- 
law potential and a subdominant energy density in a matter (radiation) dominated 



universe, Eq. (|3.8|) possesses |5(J an exact scaling solution, which behaves as a sta- 
ble attractor provided that the power is greater than 6 (10). In this case, the field 
is eventually trapped in the attractor and its energy density decreases relative to 
the dominant matter (radiation) energy density as time elapses. Moreover, any pri- 
mordial perturbations in this field are washed out since the field falls in the same 
attractor no matter what its initial conditions were. Consequently, under these cir- 
cumstances, the field could not play the role of a curvaton since the limiting would 
come out tiny. Fortunately, all the powers involved in our PQ potential in Eq. (|3.1| ) 
are less than or equal to 6 and, thus, no attractor type behavior is encountered with 
this potential as our numerical findings also clearly show. Furthermore, if our full 
potential V(0) in Eq. (|3~7| ) is dominated by the variable induced mass in Eq. ( |3.6| ), 
new scaling type solutions appear |T7|]. However, the initial perturbations are not 
washed out in this case. 



- 22 - 




Figure 5: The evolution of the field <fi for two (pfs differing by a small amount and lying in 
the 3rd band near its (a) lower, or (b) upper boundary at A = 1CP 4 . The magenta (brown) 
field paths correspond to <f)f being nearer (further from ) this boundary. The maximum at 
4> = —(praax acts as a 'tachyonic amplifier' of the difference between the fields. 

As already indicated, the limiting comes out generally small. However, if (f>f is 
near the upper or lower boundary of the 3rd band, can be enhanced considerably. 
To see this, let us consider two <p/s slightly above (below) the lower (upper) boundary 
of this band which differ by a small amount 5(f). We take again A = 10~ 4 . The relevant 
parts of the resulting field paths are shown in Fig. [5]a (Fig. |5|b). The magenta path 
corresponds to the <f>f which is nearer the boundary, while the brown one to the 0/ 
furthermost from the boundary. We observe that the 'magenta field' is stabilized 
nearer the maximum at <f> = — max than the 'brown field' after its first slow (its slow 
return from its first slow) oscillation to negative values. Then, due to the flatness of 
the potential around the maximum, the 'magenta field' remains near the maximum 
longer than the 'brown field'. Consequently, it falls towards the PQ minimum at 
= —f a with a relative delay and, thus, the gap 5(f) between the fields increases 
considerably. This leads to a sizeable enhancement of the disparity in the amplitudes 
of the subsequent damped oscillations of the two fields about the minimum. 

We see that the local maxima can act as 'tachyonic amplifiers' of the curvaton 
perturbation which originates from inflation. Note, however, that the amplification 
happens only if 0/ is close to the boundary of the band and only once during the field 
evolution, since, at subsequent times, the field never approaches the local maxima 
again. This effect is crucial for the viability of our model since otherwise the limiting 
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turns out too small. It should be mentioned that an amplification of perturbations 
while the field is near a local maximum of the potential was also encountered in 
Ref. . However, the phenomenon there was of a quantum nature referring to the 
quantum perturbations of the inflaton in de Sitter space. In our case, the effect is 
completely classical and takes place well after the termination of inflation. 

5.4 The curvaton decay 

We now turn to the discussion of the curvaton decay and the evaluation of the final 
curvature perturbation. The curvaton decays into a pair of Higgsinos via the second 
coupling in the superpotential of Eq. ( |2.2| ). This same coupling is also responsible 
for the /x-term. Of course, for this decay to be kinematically possible, we must make 
sure that the Higgsino mass /i = j3f 2 /4mp does not exceed half of the curvaton mass 
which is given by 



ml 



^' Q (0 = ±/ B ) = ^\A\ 2 -12 (\A\ + vW 1 ^) rn\ /2 . (5.6) 



This mass is A-independent and equals about 965 GeV for the input numbers used. 
The curvaton decay time = FT , where is the curvaton decay width given by 

(3 2 f 2 

T * = 8^F P m *- (5 ' 7) 

We take ll = 300 GeV, which yields (3 ~ 0.697A and ~ 7.6 A" 1 x 10 13 GeV" 1 . 

It is important to ensure that the coherently oscillating curvaton field does not 
evaporate [^j as a result of scattering with particles in the thermal bath before it 
decays into Higgsinos. During the period between the onset of field oscillations about 
a PQ vacuum and roughly the electroweak phase transition, the Higgs superfields are 
in equilibrium with the thermal bath and could knock the PQ field out of the zero 
mode. The dominant process is the Higgsino-curvaton scattering through the second 
coupling in the RHS of Eq. ( |2.2| ). The cross section is a ~ j3 2 /m 2 P , which yields the 
thermal decay rate Tth ~ /3 2 T 3 /m 2 P . Avoidance of evaporation then requires that 
(3 2 T s /m 2 P < H. Before (after) reheating, H ~ T 4 /m P T 2 eh |3§ (H ~ T 2 /m P ) and 



this requirement becomes T > /3 2 T 2 eh /m P (T < m P /f3 2 ), which is well-satisfied. 

The energy density of the oscillating curvaton at is evaluated numerically for 
each 0/ in the 3rd band and each A in the range 10~ 4 — 10 -3 . From Eq. ( |4.5| ), we can 
then calculate /dec, the fraction / at curvaton decay. The final curvature perturbation 
( is estimated, in the instantaneous decay approximation, by multiplying this fraction 
by the limiting (see Eq. (|4.6|)). The result is put equal to 5 x 10~ 5 as required 
by COBE [p9| . We find that there are two solutions for each A, one close to the 
lower boundary of the band and the other close to the upper boundary. These 
solutions lie on the two blue lines in Fig. where the green and red areas together 
constitute the 3rd band. It is true that the solutions turn out to be quite close to 
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Figure 6: The 3rd band (all the colored area) with the blue solution lines on which the 
COBE constraint ( = 5 x 1CP 5 is satisfied. The green (red) area of the band is allowed 
(excluded) by the non-Gaussianity bound from the WMAP satellite. 

the boundaries of the band. This means that the 'tachyonic amplification' of the 
inflationary perturbation in the curvaton field, which occurs near a maximum of the 
potential, plays a crucial role. On the other hand, the gap between the solutions and 
the boundaries is not so tiny as to imply that the solutions are finely tuned. Moreover, 
this gap is huge compared to the inflationary perturbation 5<ft — 3.14 x 10 10 GeV. 

We already stressed that, in the curvaton scenario, significant non-Gaussianity 
of the curvature perturbation appears if the curvaton decays before dominating the 
energy density of the universe. The recent CMBR data from the WMAP satellite 
yield |J an upper bound on the possible non-Gaussianity in the curvature pertur- 
bation. This bound implies |Tl[ that, at 95% c.L, the fraction f dec > 9.33 x 10~ 3 . 
The corresponding allowed (excluded) area of the 3rd band is colored green (red) in 
Fig. |6|. We see that the (blue) solution lines lie totally in the green area and, thus, are 
fully consistent with the negative results of the WMAP satellite on non-Gaussianity. 

As already explained, if 0/ lies in a green band and near one of its boundaries, 
the PQ field, before starting its damped oscillations about a PQ vacuum, stays for 



a while near a maximum of the potential. This 'tachyonic effect' is also important 
for the fulfillment of the non-Gaussianity constraint. Indeed, it causes a delay in 
the commencement of the damped oscillations of <fi and, thus, leads to a larger 
at curvaton decay. Therefore, the fraction /d ec is enhanced. As a consequence, the 
non-Gaussianity bound can be satisfied provided that 0/ is adequately close to the 
boundaries of the band. This explains the fact that the green area in Fig. ^ extends 
near the boundaries of the band. Needless to say that the enhancement of /dec also 
facilitates the fulfillment of the COBE constraint on the curvature perturbation. 

In summary, we conclude that, in the investigated scheme which simultaneously 
solves the strong CP and p problems, the PQ field can successfully and naturally act 
as a curvaton with all the relevant cosmological requirements satisfied. 

6. Conclusions 

We considered a minimal extension of MSSM which simultaneously solves the strong 
CP and p problems via a PQ and a continuous R symmetry. This model can be 
readily embedded in simple and natural SUSY GUTs which incorporate the hybrid 
inflationary scenario and generate the observed baryon asymmetry of the universe 
through a primordial leptogenesis. We, thus, took the model supplemented with 
hybrid inflation and leptogenesis, but without committing ourselves to the specific 
details of these scenarios. 

We examined whether, in this model, the PQ field can play the role of the curva- 
ton generating the primordial density perturbations which are needed for explaining 
the structure formation in the universe and the observed anisotropies in the CMBR. 
A crucial requirement is that the PQ field is light during the last 50 — 60 inflationary 
e-foldings so that it receives a superhorizon spectrum of perturbations from inflation. 
This condition can be easily satisfied in our model with global SUSY, where the PQ 
potential is almost flat, provided that, during inflation, the value of the PQ field is 
not too large. We must, though, employ some mechanism which can prevent the 
lifting of the flatness of the potential by the SUGRA corrections during inflation. 

We followed the evolution of the PQ field after the end of inflation under the 
assumption that the SUSY breaking corrections to its potential from the finite energy 
density of the early universe are somewhat small. We found that the field, after an 
initially slow evolution, is set in damped oscillations about zero or one of the PQ 
vacua, depending on its value right after inflation. More precisely, the values of the 
PQ field at the end of inflation can be classified into successive bands which lead to 
the trivial or the PQ vacua in turn. Needless to say that only the bands leading to 
one of the PQ vacua are useful for our purposes. 

In our model, as it turns out, the bulk of the CDM in the universe consists of 
axions which are produced at the QCD phase transition. They carry an isocurvature 
perturbation which is uncorrelated with the total curvature perturbation generated 
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by the PQ field. The rest of CDM is made of LSPs (lightest neutralinos) which 
originate from the late decay of the primordial gravitinos produced thermally at re- 
heating. These neutralinos as well as the baryons, which are assumed to come from a 
primordial leptogenesis that took place at reheating, acquire an isocurvature pertur- 
bation fully correlated with the curvature perturbation. So, the overall isocurvature 
perturbation has a mixed correlation with the total curvature perturbation. We 
found that the presently available restriction on such an isocurvature perturbation 
from the CMBR and other data excludes the lowest lying band which leads to a PQ 
vacuum. High order bands are also excluded since they do not fulfill the condition 
that the curvaton is light during inflation. We are, finally, left with only one allowed 
band leading to a PQ vacuum. We, thus, focus our attention to this particular band. 

The inflationary perturbation of the PQ field evolves as the cosmic time elapses 
after the end of inflation and, when the PQ field is set in harmonic oscillations about 
a PQ vacuum, yields a stable perturbation in the energy density of this field. At 
the decay of the PQ field, the perturbation is transferred to the radiation dominated 
plasma generating the total curvature perturbation. We found that the COBE con- 
straint on the total curvature perturbation can be naturally fulfilled provided that we 
are somewhat near the boundaries of the band. An important phenomenon, which 
allows us to achieve this, is the 'tachyonic amplification' of the original inflation- 
ary perturbation in the PQ field as this field is temporarily stabilized near a local 
maximum of the potential. 

Finally, we imposed the bound on the possible non-Gaussianity of the total 
curvature perturbation from the recent CMBR data obtained by the WMAP satellite. 
We saw that, although the bulk of the band is excluded by this constraint, the COBE 
requirement is still fully consistent with it. The 'tachyonic effect', i.e. the fact that 
the field hangs around a maximum of the potential for a period of time if its value 
right after inflation is close to the boundaries of the band, is crucial for having the 
non-Gaussianity requirement satisfied. The reason is that, due to this effect, the 
commencement of the damped field oscillations is delayed and, thus, the density 
fraction of the curvaton is enhanced. This also facilitates the achievement of the 
COBE constraint on the curvature perturbation. 

In conclusion, we have shown that, in our scheme, which naturally and simulta- 
neously solves the \x and strong CP problems, the PQ field can successfully play the 
role of the curvaton generating the total curvature perturbation in the universe in 
accord with the COBE measurements. The bounds on the isocurvature perturbation 
and non-Gaussianity from the available CMBR data can also be satisfied. 
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